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ABSTRACT 

Detailed models are compared to recent infrared observations of the nearby extrasolar planet, HD 
189733b. It is demonstrated that atmospheric water is present and that the planet's day side has 
a non-isothermal structure down to gas pressures of ~ 0.1 bars. Furthermore, model spectra with 
different amounts of CO are compared to the observations and an atmosphere absent of CO is excluded 
at roughly 2-sigma. Constraining the CO concentration beyond that is unfortunately not possible 
with the current Spitzer photometry. However, radically enhanced (or depleted) metal abundances 
are unlikely and the basic composition of this planet is probably similar to that of its host star. When 
combined with Spitzer observations, a recent ground-based upper limit for the K-b&nd day side flux 
allows one to estimate the day-to-night energy redistribution efficiency to be ~ 43%. 

Subject headings: planetary atmospheres - extrasolar planets 
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1. INTRODUCTION 

Hundreds of extrasolar planets have now been discov- 
ered along with over a dozen transiting systems. Despite 
the rapid pace of discovery, many important questions 
remain unanswered, especially concerning the transiting 
giant planets (hot-Jupiters). How efficiently is incident 
stellar energy absorbed by the day side redistributed to 
the night side? Do hot-Jupiters have atmospheres with 
enhanced metal abundances relative their host star? And 
under what conditions (and depths) are the day side 
atmospheric temperature structures isothermal? For- 
tunately, for the handful of nearby transiting planets, 
direct flux measurements have been obtained using the 
Spitzer Space Telescope. Fluxes provide a direct look into 
a planet's depth-dependent atmospheric conditions and 
perhaps hold the answers to some of these questions. 

HD 189733b , is the closest kno wn jovian-mass tran- 
siting planet (|Bouchv et al.l 120051 ). and flux measure- 
ments of its day side are available across a wide range 
of infrared wavelengths along with a re cent measure- 
ment of the planet's night s i de at 8 nva. (Domin g et al.l 
20061 iGrillmair et all 120071 : ICharbonneau et all l200a 
Knutson et al J 120071) . HD 189733b is expected to have 
photospheric temperatures in the 1000 to 1500K range 
on its day side maintained by the large amount of inci- 
dent stellar flux. At just 20 pc away from Earth, this 
planet is a valuable test case for atmospheric dynam- 
ics and chemical enrichment by accretion of extrasolar 
comets and asteroids. 

In this Letter atmospheric models are compared to re- 
cent IR observations of HD 189733b to infer some of 
its basic photospheric properties. The concentrations 
of important molecules (e.g., water and CO) are dis- 
cussed along with the day-to-night energy redistribution 
by depth-dependent horizontal circulations. 

2. ECLIPSE DEPTHS AND SPECTRA 

As a transiting planet passes behind its host star a 
small but measurable drop in infrared flux occurs as the 
planet's day side contribution disappears. A few hours 
later, the planet reappears from behind the host star and 



the total system flux returns to its previous value. The 
change in flux during this secondary eclipse provides a 
precise measurement of the day side planet-star flux den- 



sity ratio (e\ 



RiF* 



where R p .± are the planet and 



star radii) at a given wavelength. Since even heavily irra- 
diated planets are much cooler than their host star, their 
peak flux emerges at a significantly longer wavelength 
(\™ ax ) than it does for the host star causing e\ to rise 
toward longer wavelengths and making secondary eclipse 
observations more practical at A > X£ ax . Regardless of 
the detailed absorption features in the two spectra, e\ 
falls rapidly from the infrared to the optical as a result of 
the steep rise of the stellar spectrum and, consequently, 
the deepest eclipse measurements made with Spitzer are 
generally where the planet flux density is small. 

Secondary eclipses of HD 189733 have been observed 
photometricall y with Spitzer at 3 . 6, 4.5, 5.8, 8.0, 16 
and 24 mn foeming et all 120061 : ICharbonneau et all 
120081: iKnutson et alj 120071) A secondary eclipse has 
also been observed sp ectroscopically from 7.5 to 15 /an 
(jGrillmair et al.ll2007l ). Figure Q] shows all of the day 
side Spitzer secondary eclipse photometry for HD 189733 
compared to a synthetic e\ spectrum produced by di- 
viding a synthetic planet spectrum by a synthetic stel- 
lar spectrum (optimized for the K-type stellar host). 
The e\ spectrum has also been scaled by the planet- 
star surface areas d etermined from transit observations 
(|Bouchv et alJl2005t ). 

Given that e\ emphasizes the faintest portions of the 
planet's spectrum, comparing actual planet fluxes to 
model spectra is more revealing. The observed fluxes 
shown in Fig. [3] were calculated using the out-of-eclipse 
stellar fluxes and the definition of above. As can 
be seen, the planet's spectrum peaks at ~ 3.6/xm - far 
from the deepest Spitzer eclipse depth. Also, while e\ 
is nearly flat from 3.6 to 4.5 /im, the individual planet 
spectrum is quite steep across this wavelength range. 

Overall, the observations agree fairly well with the 
model shown in Fig. [2] (same model in Fig. [I}. The stan- 
dard model includes a solar abundance composition in 
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Fig. 1. — Model planet-star flux density ratios (e\) assuming 
no energy is redistributed to the night side. White circles are 
the Charbonncau ct al. ( 2008) Spitzer secondary eclipse measure- 
ments with 1-<t error- bars. White boxes are additional 8.0 an d 
16 )im eclipse values (Doming ct al. 2006; Knutson ct al. 2007). 
White diamonds are the model values integrated across the Spitzer 
instrumental response curves, indicated by dotted lines. 



chemical equilibriu m and irradia t ion fo l lowin g the meth- 
ods described in iBarman et all (|200lL 120051 ). Further- 
more, this model assumes no redistribution of absorbed 
stellar flux from the day side to the night side and, 
hence, implies a very cold night side. Interestingly, zero 
energy redistribution to the night side is at odds with 
recent 8 (im observations indicati ng a w arm night side 
photo sphere (jKnutson et al.lf2007h . Also. iBurrows et all 
(2006) found that their full-redistribution model (50% of 
the incident energy carried to the night side) is in good 
agreement with the lDeming et all (|2006[ ) 16 /Ltm Spitzer 
observation. However, the Burrows et al. model predicts 
stronger absorption across the central two IRAC bands 
than is observed. 

3. EVIDENCE FOR WATER AND CO 

For the temperatures and pressures expected in the 
atmospheres of irradiated giant planets, models predict 
that H2O should be the third or fourth most abundant 
gas-phase molecule. With its broad absorption features 
across the IR, H2O plays an important role in sculpt- 
ing a giant planet's emergent spectrum and temperature 
structure. Therefore, measuring the water abundance 
is an important step toward understanding the planet's 
overall chemistry and structure. 

Given that basic predictions of hot-Jupiters include a 
copious water supply, the first spectra of hot-Jupiters 
were somewhat surprising in that no water absorption 
was detected. These e\ spectra turned out to be fairly 
flat from 7.5 to 11 /im which was interpreted as an 
absence of water abs o rption in the planet's spectrum 
(|Grillmair et all 120071 : [Richardson et al.ll2007l) . A few 
months after the IRS spectra were published, evidence 
for water absor ption in the at mosphere of HD 209458b 
was identified (Barman 20 0?]) and a sim ilar claim fol- 
lowed for HD 189733b (jTinetti et al.ll2007l ). These water 
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Fig. 2. — Top: Synthetic spectru m of the solar abundance model 
shown in Fig. [T] Observations (Charbonncau ct al. 2008) arc 
shown with 1-<t error-bars along with a 1450K black body (dot- 
ted line) and a model with no water opacity (in blue). Bottom: 
Standard model (red) compared to models without CO (blue) and 
CO mole fraction equal to 100 times that of the standard model. 
Diamonds are the band integrated model points at 4.5 /im. 



detections relied on observations during primary transit, 
and thus relate to different regions of the planets than 
the secondary eclipse spectra (and photometry). 

Interestingly, the Spitzer IRAC data (3.6 to 8 /xm) for 
HD 189733 show evidence of H2O and, to a lesser degree, 
CO absorption at the level predicted by model atmo- 
spheres (the mole fractions of CO and H2O for the stan- 
dard solar abundance model presented here are 6 x 10 
and 5 x 10~ 4 respectively) . The IRAC colors in particular 
are consistent with strong water absorption from ~ 4 to 8 
/Ltm. The good agreement between the no-redistribution 
model and the observations (Figs. [T]andE]) supports this 
interpretation. Furthermore, the IRS spectrum for HD 
189733 reported to have no water absorption, is reason- 
ably consistent with the model everywhere ex cept at the 
blue and red edges (|Fortnev k, Marlevl l2007h . A flat e A 
spectrum across 8 /im implies a very sharp rise in planet 
flux (see Fig. [2]) that is inconsistent with the IRAC ob- 
servations. 

Spectra with different amounts of CO are compared to 
the observations in the lower panel of Fig. [2] An atmo- 
sphere absent of CO is excluded at roughly 2-a and thus 
its presence is only tentatively suggested by the data. 
Unfortunately, the current precision of the 4.5 /jm IRAC 
photometry is not sufficient to test large departures of 
Carbon from solar ab undance or departures of C O from 
chemical equilibrium (|Cooper fc Showm an 2006). 

4. AN L-TYPE PLANET 
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Fig. 3. — Top: The synthetic planet spectrum compared to a 
real IR spectrum of L5 brown dwarf 2MASS 1507. Spitzer pho- 
tometry for HD 189733b and 2MASS 1507 are also shown, scaled 
to match at 3.6 jim. Bottom: Temperatures and pressures for 
no-redistribution (top solid line) and full-redistribution (bottom 
solid line) models compared to 1250 and 1450K brown dwarf mod- 
els (dashed lines). The Spitzer and near-IR photospheric depths 
are indicated by thick lines and filled symbols. The vertical dash- 
dotted line is the approximate boundary between efficient (higher 
P) and inefficient (lower P) horizontal circulation. The dotted 
curve denotes the region of equal CO/CH4 equilibrium concentra- 
tions. 



Many hot-Jupiters have day side photospheric tem- 
peratures and pressures comparable to those in L-type 
brown dwarfs and, despite potential differences in for- 
mation and evolution, a great deal can be learned by 
comparing observations of the former to observations of 
the latter. Figure [3] compares the HD 189733b obser- 
vations and no-redistribution standard planet model to 
IRTF and Spi tzer IR observations of an L5 brown dwarf, 
2MASS 1507 (jCushing et al.ll2006l ). each scaled to match 
at 3.6 j«m. The IRAC colors of the planet and brown 
dwarf are nearly identical. The decrease in flux across the 
shortest two IRAC bands is a hallmark of CO and H2O 
absorption in brown dwarf atmospheres, lending support 
to the conclusions of the previous section. The observed 
IRS spectrum of the brown dwarf is also similar to the 
planet model spectrum, though a bit steeper. Also, the 
Spitzer photosphere in the brown dwarf is at lower pres- 
sures than in the planet where the two temperature struc- 
tures have similar slopes and produce absorption features 
with similar relative depths. 

The spectral similarities between HD 189733b and 
2MASS 1507 answers an important question - this 
planet's temperature structure can not be isothermal at 
the photospheric depths probed by Spitzer. This rules 
out one of the more plausible explanations for why the 
first set of IRS spectra showed no observed water ab- 



sorption features. Looking back at Fig. [2 one can see 
that indeed the slope of the IRS planet spectrum, more 
or less, follows that of a black body and, by itself, is 
very suggestive of a purely isothermal structure. How- 
ever, the 3.6 /jm IRAC measurement, in concert with 
the longer wavelength data, rules out any single black 
body spectrum at several a. Overall, the > 3 spec- 
tral region of the planet appears to be remarkably sim- 
ilar to that of an L5 brown dwarf believed to have a 
solar atmospheric composition dominated by water and 
a steep temperature structure. However, as discussed 
below, irradiation driven atmospheric circulations most 
likely cause the short-wavelength spectrum to depart sig- 
nificantly from that of an L5 brown dwarf. 

5. BALANCING THE BUDGET 

Assuming HD 189733b is tidally locked, the inci- 
dent stellar flux always enters the planet through the 
same hemisphere. This asymmetric heating may drive 
strong atmospheric circulations that transport significant 
amounts of energy to th e night side (|Showman fc GuillotJ 
I2002t ICho et alj|2003j) . Determining the efficiency of 
energy redistribution across the day and night sides of 
a planet is, therefore, important for un derstanding the 
global atmospheric circulation. Recently iKnutson et al.l 
(2007) measured the night side flux to be ~ 64 % of the 
day side flux at 8 /Ltm. This suggests that some fraction 
of the absorbed incident flux is, in fact, transported to 
the night side. How then can the day side Spitzer obser- 
vations so closely match the predicted flux from a model 
with zero horizontal energy transport, as shown above? 

After a few 100 Myrs, internal energy leftover from 
the formation of HD 189733b has mostly been lost, leav- 
ing the intrinsic luminosity ~ 20, 000 times smaller than 
what the planet receives from the star. Thus, for old 
short-period planets like HD 189733b, the star is the pri- 
mary energy source. Energy conservation requires that 
the planet's atmosphere radiate as much energy as it re- 
ceives and, if HD 189733b has a small albedo, the to- 
tal luminosity one should expect from the planet is de- 
termined entirely by the stellar and orbital properties; 
h to t ~ 2.6 x 10 -5 Lq. Consequently, the true bolometric 
luminosity of the planet's day or night side (Lday or 
Lnight) is needed to precisely determine the efficiency of 
redistribution. Note that the efficiency of redistribution 
is difficult to constrain with single band-pass photome- 
try, since little leverage on the bolometric luminosity is 
obtained, even with a complete phase curve. 

Fig. [5] shows vF v spectra for the full and no- 
redistribution planet models. Integrating the no- 
redistribution model spectrum (top curve) for A > 3.2 
/zm accounts for 33% of L tot . Constraining \jd ay further 
requires measurements at shorter wavelengths. Fortu- 
nately, recent ground based observations have provided 
an upper limit for the day side if-band flux (less than 
1.5 mJy) that helps com plete the energy balance pic- 
ture ([Barnes et al.ll2007h . Though \™ ax ~ 3.6 /im, the 
peak of the energy distribution sits at if-band mak- 
ing the i-T-band flux a very useful gauge of the short- 
wavelength fraction of l^day By combining the day side 
measurements, one can safely conclude that a sizable 
fraction of the energy (> 30%) is redistributed to the 
night side. Furthermore, these data tell us that a 1-D 
no-redistribution model over estimates the flux at short 
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Fig. 4. — Top: planet energy distributions assuming full- 
redistribution (lower red line) and no-redistribution (top blue line) 
along with a 1450K blackbody spectrum. Over plotted are iR 
ground and space-based flux measurements for HD 189733b (see 
text for references). T he lower point at 8 /i m is the night side 
flux measurement from Knutson ct al. (2007). Bottom: pressure 
at optical depth unity for both planet models. The location of the 
Spitzer and near-IR band-integrated pressures are labeled. 



wavelengths (< 2.5 /im). 

The estimate of en ergy redistributed can be further 
improved by using the lKnutson et all (|2007l ) 8 /im night 
side measurement and the fact that only ~ 66% of L tot re- 
mains to be shared across the entire night side spectrum 
and the day side spectrum not covered by Spitzer. The 
observed night side flux is close to the full-redistribution 
model spectrum shown in Fig. 0] suggesting that the 
night side Spitzer photospheric temperatures are close 
to those predicted by the full-redistribution model. In- 
tegrating this spectrum (slightly scaled down to match 
the 8 /im night side observation) across the complete 
wavelength range accounts for ~ 43% of L tot , leaving 
~ 24% of L^t for the day side spectrum at A < 2.5 /im. 
The day side if-band upper-limit is also close to the full- 
redistribution model spectrum suggesting that both day 
and night side near-IR spectra have similar flux levels. 

As already discussed, the day side structure can not 
be entirely isothermal. The Spitzer data probe ap- 
proximately the 0.1 to 0.01 bar region on the day side 
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while </, H, and K-bands probe higher pressures, > 
100 mbars (see Fig. HJ. Lower fluxes on the day side 
at short wavelengths (as indicated at A'-band) com- 
pared to longer wavelengths suggests that energy re- 
distribution is highly depth-dependent and occurs fairly 
deep in the atmosphere (> 0.1 bars). Similar depth- 
dependent redistribution has been explored for HD 
209458b and is the expected behavior in many hot - 
Jupiters ()Cooper fc Showman! [20051 : ISeager et al.ll2005f ). 
Thus, the day and night sides probably have similar tem- 
peratures at P > 0.1 bars which are close to the full- 
redistribution structure shown in Fig. [3] The steady 
supply of incident flux from the star helps maintain the 
hot temperatures in the upper layers of the day side while 
on the night side temperatures at P < 0.1 bars are prob- 
ably more similar to those of a cooler (T e ff ~ 1250K) 
brown dwarf structure (see Fig. [3]). With lower temper- 
atures, CH4 begins to compete with CO as the dominant 
carbon bearing molecule on the night side, especially in 
the upper atmosphere. An important prediction of this 
scenario is that one should not expect significant phase 
variations at J, H, and if-band for this planet. Also, 
based on the models described here, the 24 pm night 
side flux is predicted to be ~ 75 % of the day side flux, 
a prediction that will soon be tested by Spitzer observa- 
tions. A more rigorous constraint on redistribution will 
require additional observations, especially of the night 
side, along with detailed global circulation models. 

6. CONCLUSIONS 

HD 189733b has an atmosphere rich in water and, most 
likely, CO. This planet's atmosphere also has day side 
temperatures, pressures, and chemical composition sim- 
ilar to a typical mid-L type brown dwarf. The fraction 
of energy redistributed from the day to night side has 
been estimated to be ~ 43% of the incident stellar flux. 
Furthermore, the bulk of the energy transport likely oc- 
curs at pressures higher than ~ 0.1 bars, indicating that 
the energy redistribution mechanism is most efficient at 
high pressures. It has also been demonstrated that the 
photospheric depths probed by Spitzer do not reach the 
isothermal layers predicted by atmosphere models. Con- 
sequently, an isothermal structure can not explain the 
flat e\ spectrum recently measured by IRS. 
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